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We show that it is possible to use the spatial quantum correlations present in twin beams to 
extract information about the shape of a mask in the path of one of the beams. The scheme, based 
on noise measurements through homodyne detection, is useful in the regime where the number of 
photons is low enough that direct detection with a photodiode is difficult but high enough that 
photon counting is not an option. We find that under some conditions the use of quantum states 
of light leads to an enhancement of the sensitivity in the estimation of the shape of the mask over 
what can be achieved with a classical state with equivalent properties (mean photon flux and noise 
properties). In addition, we show that the level of enhancement that is obtained is a result of the 
quantum correlations and cannot be explained with only classical correlations. 



I. INTRODUCTION 



II. BASIC CONFIGURATION 
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The field of quantum imaging studies the spatial distri- 
bution of the quantum correlations in optical fields [l|, y] ■ 
It offers the possibility of improving the resolution of 
an imaging system [3j, enhancing image detection [J] or 
beam positioning [5[ , and implementing parallel quantum 
information encoding [6| and processing [7|. 

An important goal of quantum imaging is to use the 
multi-spatial-mode properties of quantum fields to ex- 
tend the enhancement in sensitivity that can be obtained 
with quantum states from the temporal [8| to the spatial 
domain. In this paper we propose a scheme that achieves 
this goal. In particular, we show that it is possible to take 
advantage of spatial quantum correlations between two 
entangled beams to extract information about the shape 
of a mask in the path of one of the beams. We con- 
sider a configuration based on noise measurements and 
homodyne detection. As we are not extracting informa- 
tion from the mean field values, this technique becomes 
useful when the number of photons that is used to inter- 
rogate the mask is in the regime where it is low enough 
that it cannot be readily detected with a photodiode, but 
high enough that photon counting is not an option. Even 
though in this regime the mask is interrogated by a small 
number of photons, the use of homodyne detection makes 
it possible to perform the measurement. 

We analyze whether or not there is an advantage in the 
use of quantum states over classical ones with equivalent 
properties (mean photon flux and single beam noise), and 
find that under some conditions an enhancement of the 
sensitivity in the estimation of the shape of the mask can 
be obtained with the quantum states. We show that the 
enhancement that is obtained is a result of the quantum 
correlations and cannot be explained with only classical 
correlations. Finally, we find that the level of enhance- 
ment increases with the number of spatial modes in the 
field probing the mask. 
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The possibility of extracting spatial information about 
a mask from noise measurements of a beam that passes 
through it results from the fact that in homodyne detec- 
tion the shape of the local oscillator (LO) effectively acts 
as a spatial filter that selects out the portion of the field 
to be measured [9|, [l(|. That is, the spatial overlap be- 
tween the beam that goes through the mask and the LO 
affects the efficiency of the measurement. Since the beam 
acquires the shape of the mask after going through it, we 
can change the shape of the LO when performing the ho- 
modyne detection to optimize the measured noise level 
and thus estimate the overlap. Once the optimization 
criterion is met, the shape of the mask can be deduced 
from the shape of the LO. In practice, as a result of the 
type of measurement being done, this problem reduces to 
one of parameter estimation. For this particular problem 
the quantity to estimate is the overlap or mode match- 
ing between the LO and the beam that goes through the 
mask. 

The basic configuration for such a scheme is shown in 
Fig. [TJ Entangled twin beams are used to probe the mask 
whose shape we want to estimate. One of the beams goes 
through the mask and is then detected with a homodyne 
detector while the other one is directly sent to a second 
homodyne detector. The noise in the difference signal 
of the detectors is used to estimate the overlap and thus 
extract information about the shape of the mask. When 
the entangled beams that are used contain multiple spa- 
tial modes they also contain spatial quantum correlations 
which are characterized by a coherence area [ll|, UM that 
sets the spatial extent of the smallest correlated areas 
between the two beams. If the beam contains multiple 
coherence areas, and the size of these coherence areas is 
smaller than the size of the mask, then the sensitivity 
in the estimation can be enhanced. Since a coherence 
area in one beam is only correlated with a single corre- 
sponding coherence area in the other beam [13l . Il4j . it 
is important that the LOs of the two homodyne detec- 
tors select out or measure the corresponding correlated 



regions in the two beams. 
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and that this enhancement cannot be explained assum- 
ing only classical correlations. 

We compare the two-beam strategy presented in Fig.Q] 
with two analogous classically obtainable configurations. 
In the first one we consider a single beam scheme in which 
the beam that goes through the mask shares the same 
photon flux and mode structure as the entangled beams 
in the two-beam strategy. After going through the mask 
the beam is sent to a homodyne detector. The result- 
ing noise measurements are then used to estimate the 
shape of the mask. Throughout the paper we refer to 
this scheme as the single-beam strategy. In the second 
configuration we duplicate the two-beam strategy using 
classically correlated beams also sharing the same photon 
flux and mode structure as above. 



FIG. 1. Configuration for estimating the shape of a mask through 
noise measurements of spatially multi-mode entangled beams. LO: 
local oscillator; BS: 50/50 beam splitter; HD: homodyne detector; 
SA: spectrum analyzer. 



In order to implement this scheme one would like to 
perform a series of noise measurements with a dynami- 
cally changing LO. This could be done by using spatial 
light modulators to control the spatial profile of the LOs 
for beams a and b. The shape of the LOs would be mod- 
ified until the noise level is optimized, at which point 
the overlap between the LO of beam a and the field that 
went through the mask is also optimized. This implies 
that the shape of the LO of beam a corresponds to the 
one of the field that went through the mask and makes 
it possible to infer the shape of the mask from the final 
shape of this LO. 

This strategy is reminiscent of ghost imaging [15| |. in 
which the goal is to image an object without it being 
in the path of an imaging camera or detector. In ghost 
imaging the two beams that are used are composed of 
streams of correlated photon pairs. The object is placed 
in the path of one of the beams, which is then detected by 
a non-imaging, or "bucket" , detector, which alone reveals 
no information about the shape of the object. The beam 
that does not go through the mask is sent directly to a 
charge-coupled device (CCD) camera or scanning point 
detector. Correlations between the arrival of photons on 
the detectors can then be used to reconstruct the object. 
Even though originally it was believed that spatial quan- 
tum correlations were required for ghost imaging [16I418J , 
it was later realized that classical correlations could also 
be used for the same purpose [13, [2(| . As we will show 
here, there are two main differences between the scheme 
proposed here and ghost imaging. The first is that, with 
the proposed configuration, information about the ob- 
ject can be obtained directly from measurements of only 
the beam that goes through the mask. The second, and 
most important one, is that the use of quantum corre- 
lations between the beams can lead to an enhancement 
of the sensitivity in the estimation of the shape of the 
object over a classical state with equivalent properties, 



III. MULTI-SPATIAL-MODE PROPERTIES 

In order to analyze both the two-beam and single-beam 
strategies described above we need to take into account 
the multi-spatial-mode nature of the fields, even for the 
case in which only one of the spatial modes is not in 
a vacuum state. This can be done by expanding the 
field operators in terms of a complete set of spatial mode 
functions [l| , such that 



i 

b(x) =J2Mx%, 



(1) 
(2) 



where Ai(x) is an orthonormal and complete basis, x 
refers to the transverse variables of the field, and at and 
bi are the destruction operators for spatial mode i for 
beams a and b, respectively. We also assume that the 
state of the different spatial modes within each beam are 
not quantum correlated, such that the state of the field 
can be written as [21| 
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where the first N spatial modes of the field are excited 
with a wave function \ipi) for spatial mode i and the rest 
of the modes are in a vacuum state. When more than one 
spatial mode is excited (N > 1) the field is considered a 
multi-spatial-mode field. 



IV. SHAPE ESTIMATION THROUGH NOISE 
MEASUREMENTS 

The goal of the scheme described here is to estimate 
the field transfer function of the mask, which is taken to 
be of the form 



t(x) = V^We 2 



(4) 



where T(x) is the intensity transfer function and 4>{x) is 
the phase that a field would acquire after going through 



the mask and allows for the possibility of a phase mask. 
We can take into account the spatially dependent loss 
introduced by the mask by treating its effect on the field 
operator in the usual way. That is, the field operators 
are transformed according to 



a ou t(x) 



— p i<t>(x) 



y / T(x)a in (x) + \f\ - T{x)a v (x) (5) 

(6) 



bout(x) — bi n [x), 

where a v (x) is the corresponding operator for the vacuum 
field that results from losses. Note that in writing these 
expressions we have assumed that beam a is the one that 
goes through the mask while beam b has no losses. 

As described above, we are interested in obtaining 
information about the transfer function of the mask 
through noise measurements of the signal that is obtained 
with balanced homodyne detection. As is typically done 
with this type of detection, we assume the LO to be a 
classical field, such that the measured signal after beam 
a goes through the mask is given by 



In 



dx al ut (x)a LO (x) + a out (x)a* LO (x) 



= I dx 



e-^^/TW)a LO (x)al(x) 



+e-*+&y/i - T{x)a LO (x)at(x) + h.c.J , (7) 

where the integral is over the transverse profile of the 
fields, h.c. indicates the Hermitian conjugate, and 
a£o(af) is the spatial field profile of the LO. Similarly, 
for beam b the measured signal is of the form given in 
Eq. with T{x) = 1 and 0(f) = 0. 

In order to simplify Eq. (|7|) we use the mode expan- 
sion for the operators, Eq. ([1]), and expand the remaining 
product of functions in terms of the same basis, that is 



-i<f>(x) 



y/f(^jaLo(x) = aJ^Ue^Mx), 



e~^ s Wl-T(x)a LO (x) = a^t'^A^S), (8) 

i 

where a gives the amplitude of the LO for beam a and 
tie l ^ i and t' i e l< ^ i are the complex expansion coefficients 
for the product of functions in the selected basis. In par- 
ticular, ti and 4>i are the field overlap and relative phase 
between the LO, had it gone through the mask, and ba- 
sis mode i. Note that ti contains information about the 
shape of the LO and the amplitude transfer function of 
the mask, while <pi contains information about the phase 
difference between the phase of the LO and any phase 
introduced by the mask. These expansions together with 
the orthonormality property of the basis allow us to sim- 
plify I a to the form 



I a = a 



E 



uxi 



t'xl 



(9) 



where we have defined the generalized field quadrature 
for beam a before the mask, Xf 



e**oL 



-i<t>} 



The 



subindex v indicates that for those modes the field is in 
the vacuum state. 

In the same way we can simplify the expression for 
the measured signal for beam b by expanding the LO for 
beam b in terms of the basis modes Aj(x), that is 



a L o(x) = py"/ 



*Ai(x), 



(10) 



where we have assumed the LO for beam b has the same 
amplitude a as the one for beam a. In this expression at 
and 9i are the field overlap and relative phase between 
the LO and basis mode i. Note that we have defined the 
amplitude of the LO such that J2i a l — 1- We can now 
simplify if, to the form 



lb = a/ J otiX, 



b,h 



(11) 



where we have introduced the generalized field quadra- 
ture for the different spatial modes of beam b. 



A. Two-Beam Strategy 

We start by analyzing the two-beam strategy presented 
in Fig. [T] In this case entangled twin beams are used to 
estimate the shape of the mask. For this strategy, the 
signal to be measured is the noise in the difference signal, 
which is given by 



M TB = (A(I a 



h) )n, 



(12) 



where the N subindex indicates that the variance or noise 
is normalized to the corresponding standard quantum 
limit (SQL) and the TB subindex indicates the two-beam 
strategy. 

Given that beam a goes through the mask while beam 
b sees no losses, we find with the help of Eqs. ^ and (fTTj) 
that 



I a - h = a 



T.[ 



t,X. 



aiX. 



b.i 



t'xl 



(13) 



where we have assumed that the LOs for both beams a 
and b have the same amplitude a but not the same phase. 
Since we are only interested in the variance or noise of 
the signal, we can assume without loss of generality that 
the mean values for the quadratures of the field are zero, 
such that 



(/„ - h) = 0. 



(14) 



In addition, the situation when the mean value is zero 
is where the scheme described here becomes more use- 
ful, as information about the mask can only be obtained 
through noise measurements. With these expressions and 
the assumption that the different spatial modes in each 



beam are not quantum correlated, we can show that the 
normalized noise in the difference signal is given by 



Mtb — 



-21^X^X1) -t] 



(15) 



where we have used the relation J^. tf = 1 — ^L if and 
have normalized the noise to the SQL for beams a and 6, 
that is, we divide by 2|ct| 2 . 



B. Single-Beam Strategy 

In order to consider the configuration in which only 
classical resources are used, we analyze the single-beam 
strategy. In this case only the beam that goes through 
the mask is detected and the signal to be measured is the 
normalized noise of that beam. The measured signal for 
this strategy can be obtained from the one of the two- 
beam strategy, Eq. (fT5j) . by setting a, = and taking 
into account that for the single-beam strategy the nor- 
malization factor (the SQL) is given by \a\ 2 instead of 
2\a\ 2 . We thus have that 



M SB = 



= (Ap a ) N = J2tl[(Mxi 



) 2 >-l 



1. 



(16) 



where we have again assumed (I a ) = 0. 

It is important to note a coherent state cannot be 
used with this strategy to extract information about 
the transfer function of the mask, as the summation 
in Eq. (|16l) vanishes since a coherent state always has 
(A(Xf i ) 2 ) = 1. As a result, Msb = 1 and no informa- 
tion can be obtained about the coefficients ti and thus 
about the shape of the mask. This is also the case for 
the two-beam strategy. The reason that a coherent state 
is not useful is because it has the same noise level as a 
vacuum state. As a result, the vacuum noise contribution 
that is introduced by the losses due to the mask will not 
change the noise level of the field after it goes through the 
mask. Thus, this strategy requires a state whose noise 
level is initially different from that of a coherent state. 
This is an important point since for this problem a coher- 
ent state does not provide a classical limit, as is usually 
the case with other estimation techniques. 



C. Shape Estimation 

As described above, the problem being considered re- 
duces to one of parameter estimation. In particular, the 
parameter to estimate is the overlap between the LO and 
the beam that goes through the mask. We can param- 
eterize the noise of the signal that is measured with ho- 
modyne detection in terms of the transmission of the LO 
that one would obtain if it had gone through the mask 



instead, and make it the parameter to estimate. The 
transmission of the LO through the mask is then given 
by 



T = 



J dxT(x)\a LO (x)\' 
J dx\a LO (x)\ 2 



Y,tl 



(17) 



where the integral is over the transverse profile of the LO. 
In obtaining this result we have used Eq. © to express 
T in terms of U. 

For a general mask transfer function for which different 
transmission levels are possible, one would have to im- 
plement a different protocol in which a series of measure- 
ments is performed. For example, one could use a series 
of LOs with different shapes, such as the basis modes, 
in order to extract enough information to estimate the 
transfer function of the mask. This would make it possi- 
ble to extract all the values of £;. Once all these values 
are known and the input state is completely characterized 
it would be possible to reconstruct the complete transfer 
function of the mask. 

For the case of a phase mask, information about the 
phase of the transfer function is contained in the phase 
factor of the expansion coefficients of Eq. ©, that is 
fa. As can be seen from Eqs. (|T5|) and (fTl)l) this phase 
factor appears in the generalized quadrature and selects 
the quadrature that is measured. For a phase sensitive 
state in which the noise properties depend on the phase 
being measured, one can modify the phase front of the LO 
for beam a in such a way as to minimize or maximize the 
noise. Once this is done the phase of the transfer function 
of the mask can be extracted from the phase front of 
the LO. For the two-beam strategy, which requires a LO 
for each beam, independent control of the phase front 
of each LO would be required, as the LO for beam a 
will have to compensate for the additional phase that the 
beam acquires after going through the mask. It should 
be noted that a phase insensitive state cannot be used to 
obtain information about a phase mask, as the measured 
noise will be independent of fa and thus will not provide 
information about the phase of the transfer function. 

For the protocol as described here to work we require 
the intensity transfer function of the mask at each posi- 
tion x to have a transmission of either zero or one. We 
assume that the area of the LO is the same as the one of 
the mask for T — 1. We also assume that the noise sig- 
nal that is used for the estimation, Mtb or Msb, either 
grows or decreases as the transmission T increases, such 
that a maximum or minimum is obtained when T = 1. 
Under these conditions once the noise level is optimized 
we know that the LO and the field that went through 
the mask have the same shape. Thus, one can then ob- 
tain the intensity transfer function of the mask from the 
shape of the LO. 

As can be seen from Eqs. (fT5"|) and (fT6|) only spatial 
modes with a noise level different than the one of a co- 
herent state contribute to the measurement. As a result, 
even if only one spatial mode is not in a vacuum state, it 



is possible to obtain information about the transfer func- 
tion of the mask. This means that a multi-spatial-mode 
field is not required for the estimation. A single-mode 
field, however, is in general not detected with unit effi- 
ciency (t 2 ^ 1) when the LO has the same shape as the 
mask and T = 1. For an arbitrary mask, when the LO 
matches its shape, the expansion of the LO in terms of 
the basis modes requires more than a single mode. As we 
will show later, this reduction in efficiency will tend to 
lead to a reduction in the accuracy or sensitivity of the 
estimation of T and having multiple spatial modes will 
make it possible to obtain a more accurate estimate of 
T. 



SENSITIVITY OF THE ESTIMATION 



We thus find that the sensitivity in the estimation of 
T with the two-beam strategy is given by 



AT 2 B = 



2M TB 



\dM TB /dTY 



8M TB 



aE i (t?(A(jy» i ) 2 )+a?<A(x» i )='>-2t i a i <.y* i .y» i »] 

dT 

(22) 

where Mtb is given by Eq. (JTSJ). In general this expres- 
sion needs to be evaluated numerically and its behavior 
depends on the mode structure of the entangled beams 
and how the beam that goes through the mask is atten- 
uated by it. 



It is important to study the sensitivity that can be 
achieved in the estimation of the transmission T. In gen- 
eral, the uncertainty in estimating parameter P through 
measurements of quantity M is given by 



AP 2 



AM 2 



\dM/dP\ 2 



(18) 



which takes into account the slope of the measured sig- 
nal with respect to the parameter being estimated and 
the noise on the measurement. As this uncertainty is re- 
duced, the sensitivity of the estimation is increased, thus 
making it possible to obtain a better estimate for P. We 
now calculate the sensitivity for both strategies described 
above in the estimation of T. 



A. Two-Beam Strategy 

To determine the sensitivity of this strategy we need 
to calculate the noise of Mtb , which means that we need 
the variance of the difference noise signal, that is 



AM TB = ((I a - 7 6 ) 4 ) - ((/„ - i b ff 



(19) 



To simplify this expression we assume that the fields 
that are used have Gaussian statistics, which means that 
higher order moments can be expressed in terms of first 
and second order moments [22| . In particular, this im- 
plies that for the quadrature of each mode 



{{x?f) = H(x?yy 



(20) 



Using this property we find that Eq. (Tl9l) can be simpli- 
fied to 

AM TB = Ul-T + Yl[ti( A ( jt ii) 2 ) 



+ai(A(XLr)-2tMKiKi) 



2M TB . 




B. Single-Beam Strategy 

To calculate the sensitivity of the single-beam strategy 
we need to calculate the variance of M$b , which is given 
by 



AM 



SB 



<4 4 > - (ilf 



(23) 



Using the mode expansion given in Eq. © for I a and the 
Gaussian property of the state given in Eq. (|20[) we find 
that 



AM 2 SB 



£*?< A (*t 



1-T 



2M 2 B . 



(24) 



We thus have that the sensitivity for the single-beam 
strategy is given by 



ATj B 



2M| B 



\dM SB /dT\' 



= 2 



Es«?<A(je*) a > + i-r 



a(X;,^(A(x^) 2 



OT 



(25) 

We again find that an analytical solution as a function of 
T is not possible in general. Note that this result is also 
valid for quantum states such as single-mode squeezed 
states. Even though using this type of state might lead 
to an enhanced sensitivity when compared to the same 
strategy using classical states, we will only consider clas- 
sical states with the single-mode strategy in the rest of 
the paper. 



Enhancement in Sensitivity with Quantum 

States 



The uncertainty in the estimation of T given in 
Eqs. (|22"T) and ([25]) are general results that apply to any 
state with the only restriction being that it have Gaus- 
sian statistics. In general whether or not an enhancement 



in sensitivity can be obtained with the two-beam strat- 
egy using quantum states over the single-beam strategy 
using classical states depends on the actual mode struc- 
ture and shape of the mask. The functional form of these 
results will change for each configuration and, in general, 
a numerical analysis would be required to compare the 
two strategies. 

In order to compare the two strategies we consider the 
limit in which all the spatial modes have the same noise 
properties, that is, <A(A^) 2 ) = (A{Xff), <A(A^) 2 ) = 

<A(X«) 2 ), and (XtA) = {X+X°). This situation 
would correspond to the limit in which the state is com- 
posed of an infinite number of spatial modes all with the 
same properties. Even though this is not a realistic limit, 
it is still instructive to consider it, as it allows us to ob- 
tain an analytical solution for Eqs. (PI21 and (P25l) and 
illustrate the enhancement that can be obtained through 
the use of quantum states. 

For the two-beam strategy we assume that the quan- 
tum state of light used is given by vacuum twin beams. 
For this type of state each beam by itself is phase insensi- 
tive, so we can drop the phase dependence ((A(X$ ) 2 ) = 
(AX 2 ) and (A(X[) 2 ) = (AX 6 2 )). In addition, the noise 
properties of both individual beams are the same, that is 
(AX b 2 ) = (AX 2 ). In this case Eq. (JT5J) reduces to 



a state under the conditions described above we find that 



M TB = 1 + - 



(A(X a ) 2 )-l (T+l) 



<A(*J 2 >--MoH5>«*> ( 26 ) 



AT| fl = 2 



(r(AX 2 ) + 1 - Tf 
((AX 2 ) - l) 2 



(27) 



Figure [5] compares the uncertainty in the estimation 
of T for both strategies in the case in which the phases 
of the LOs are set to measure the minimum noise in the 
quadrature difference with a value of Mq = 0.1. For 
a minimum uncertainty state this implies a normalized 
single beam noise (AX 2 ) = 5. As we can see from this 



< 



Thermal state 




where we have written (XfX®) in terms of the noise 
properties of the input beams (before the mask) and de- 
fined the input normalized noise of the quadrature differ- 
ence M$' 6 = {A(X% - Xlf) N . For the correct choice of 
4> and 9, Mq corresponds to the amount of quadrature 
difference squeezing in the twin beams J23| . In order to 
simplify this result further we assume that the shape of 
the LO for beam b is not changed and that it selects out 
only one spatial mode, that is, a± — 1 and a^i = 0. 
Note that this will not lead to optimum results, as the 
LO for beam b will not be matched to the one for beam a. 
This will lead to the detection of uncorrelated coherence 
areas between beams a and b. Finally, we assume that for 
the spatial mode of beam a corresponding to the one that 
was selected by the LO in beam b the field transmission 
through the mask scales linearly with the transmission. 
As an example we take t\ — Q.ST. 

For the single-beam strategy we consider a classical 
state with the same single beam noise and mean photon 
flux as the beam that goes through the mask in the two- 
beam configuration. Such a state corresponds to a ther- 
mal state and is exactly the state that is obtained from 
the vacuum twin beams when only one of the beams is 
measured [22|]. Even though this state is obtained from 
an entangled state, it can be produced with a classical 
source and can be considered a classical state. For such 



FIG. 2. Uncertainty in the estimation of the transmission 
(AT 2 ) as a function of the transmission T for a thermal state 
(single-beam strategy) and twin beams (two-beam strategy) 
for (AXl) = 5 and M =0.1. 

figure, the two-beam strategy with twin beams offers an 
enhancement in sensitivity over the single-beam one with 
a thermal state for larger values of T (T > 0.65 for the 
particular example shown in Fig. [5]) . The value of T at 
which such an enhancement occurs depends on the initial 
noise properties of the states and the functional form 
assumed for t\. The largest enhancement is obtained for 
T = 1 which is where we want to operate, since it is 
at this value of T where the LO and the field that went 
through the mask have the same shape and an estimation 
of the shape of the mask can be made. 

We now consider the more favorable case for the two- 
beam strategy in which the LOs for beams a and b only 
measure correlated coherence areas. If we assume that 
the mode structure for beams a and b is the same and that 
quantum correlations are only present between modes 
with the same spatial profile, then this case implies that 
the LO for beam b has the same shape as the LO for 
beam a after the mask. This will only be the case for 
values of T s=a 1. For this region of operation, having 
the same shape for both LOs implies that a% ~ ti, which 



allows us to simplify Eq. (I22j) to the form 



ATj, B w 2 



(l-T)((AX 2 ) + l) + 2TA/ 



2M - ((AX 2 ) + 1) 



(28) 



In this case, the sensitivity for the single-beam strategy 
is still given by Eq. (J27I) . In the limit of large levels of 
squeezing ((AX 2 ) > 1 > M ) and T = 1 wc find that 
the sensitivity is enhanced by a factor 



ATJ B 




AT 2 



(29) 



which shows that the sensitivity enhancement factor 
grows with the level of entanglement or squeezing present 
in the twin beams. 



D. Enhancement due to the Multi-Spatial-Mode 
Nature of the Field 




FIG. 3. Enhancement in sensitivity of the estimation of T of the 
two-beam strategy with respect to the single-beam as a function 
of the number of excited spatial modes N . An enhancement is 



obtained for AT| B 



/AT| 



> 1. 



As we have shown, a multi-spatial-modc field is not 
required to estimate the transfer function of the mask, 
as a field with only a single excited spatial mode can 
provide enough information for the estimation. Next, 
we study whether there is any advantage to the use of 
a multi-spatial-mode field by reducing the problem to a 
one dimensional one and calculating the uncertainty in 
the estimation of T numerically. To do so we consider 
a square aperture with a transmission of one inside the 
square and zero everywhere else. We assume that the 
LOs for both beams have the same shape and size as the 
mask with a constant amplitude, such that the goal is to 
find the position of the mask. The problem is reduced to 
a one dimensional one by limiting the movement of the 
LO in only one direction. 

To perform the numerical analysis we assume a 
Hermite- Gauss basis for the mode decomposition and 
calculate the enhancement in sensitivity obtained with 
the two-beam strategy over the single-beam strategy, 
AT| B /AT^ B , as a function of the number of excited spa- 
tial modes N in the field. We assume that all the modes 
that are excited have the same noise properties as before, 
(AX 2 ) = (AX 2 ) = 5 and M = 0.1. Figure [3] shows the 
enhancement in sensitivity as a function of N for T = 1 . 
This value of the transmission corresponds to the one at 
which one would operate to estimate the shape of the 
mask. In this figure a value of ATg B /ATj. B > 1 indi- 
cates an enhancement in the estimation of T with the 
two-beam strategy based on a quantum state over the 
one-beam strategy based on a classical state. Due to the 
symmetry of the basis modes and the aperture, the odd 
Hermite- Gauss modes do not contribute to the estima- 
tion of T. As a result, only even terms are taken into 
account in the numerical analysis. In obtaining this re- 
sult we have added the modes from lower to higher order. 
That is, N = 1 contains only mode TEM 00 , N = 2 adds 



mode TEM 2, N — 3 adds mode TEM 2 o, and so on un- 
til all the even modes between TEMoo and TEMgs have 
been added for N = 25. For the particular mode struc- 
ture and mask we have used we find that for one and 
two spatial modes the classical strategy actually does 
better than the quantum one. For a larger number of 
modes, however, we do get an enhancement with quan- 
tum states. Additionally, we find that as the number of 
modes increases the level of enhancement in sensitivity 
also increases and saturates at an enhancement factor 
of around 3.2. This result shows that as the number of 
spatial modes present in the field that probes the mask 
increases so does the enhancement that can be obtained 
with quantum states over classical ones. 



E. Role of Quantum Correlations 

A natural question to ask is whether the enhancement 
that can be obtained with entangled states is due to the 
quantum correlations or if it can also be obtained through 
classical correlations. To answer this question we need to 
consider the two-beam strategy in the presence of only 
classical correlations. That is, we assume that the state 
of the field is described by the density matrix [24j 



P 



^PtPa,i ® Pb,; 



(30) 



where p a ,i and pb,i are the density matrices for beams a 
and b, respectively, and pi allows for a statistical mixture 
of separable states. This expression makes it possible to 
consider classical correlations. 

For the state of the field given by Eq. (|5U)l we have 



8 



that 



(*liKi) = Y,Vi&ii)i$h)i = °- ( 31 ) 



With this result, the sensitivity in estimating T given by 
Eq. (US]) reduces to 



AT^ Q 



E 4 (^(A(^) 2 )+a 2 (AX« 



1 -T 



dT 



(32) 

where the subindex NQ indicates that there are no quan- 
tum correlations between beams a and b and we have as- 
sumed (Xi to be independent of T. Since (A(X£ J 2 ) > 



we can see by comparing Eq. (1521 with Eq. (j2"5)) that 
ATjyg > ATg B . That is, the presence of only classical 
correlations in the two-beam strategy does not make it 
possible to obtain an enhancement in sensitivity over the 
classical single-beam strategy and in fact, in general, it 
leads to a decrease in sensitivity. This result shows that 
quantum correlations are needed in order to obtain the 
enhancement in sensitivity. 



VI. CONCLUSION 



the path of an optical beam. For this purpose we have 
analyzed two different strategies: a two-beam configura- 
tion based on quantum states and a one-beam configu- 
ration based on classical states. We have shown that an 
enhancement can be obtained with the use of quantum 
states over classical ones. This enhancement is a result of 
the quantum correlations present between the two beams 
and cannot be explained through the presence of only 
classical correlations. 

Even though a multi-spatial-mode field is not required 
for the estimation of the transfer function in either the 
two-beam or single-beam configurations, we have shown 
that the presence of multiple spatial modes can increase 
the enhancement in sensitivity that is obtained with the 
quantum states. In addition, the level of enhancement 
increases with the level of entanglement or squeezing be- 
tween the twin beams used for the estimation. 

Although we have limited the analysis to the case of 
a mask with a transmission value at a given point of ei- 
ther zero or one, we have shown that it is possible to ex- 
tend the analysis presented to obtain information about 
a phase mask and to treat the case of soft apertures. This 
makes it possible to estimate the complete transfer func- 
tion of an arbitrary mask through noise measurements. 



We have shown that it is possible to use noise mea- 
surements to estimate the transfer function of a mask in 
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